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Abstract — In this paper, the authors have discussed and shown 
how to tune the PID controller in closed loop with time-delay 
for the double integrator systems for a particular stability 
margins. In math model it is assumed that time delay (6) of 
the plant is known. As a case study the authors have consid- 
ered the mathematical model of the real-time beam and ball 
system and analyzed the simulation and real time response. 

Index Terms — double integrator, PID, stability, time delay 

I. Introduction 

JW atkins [1] worked with the PD control for double 
integrator systems with time delay. This paper is an extension 
in which PID control is analyzed in simulation as well as in 
real time. Integral part of the controller eliminates steady- 
state error, which can be necessary in this kind of systems. 

Equations delivered in this paper and m-files based on 
them can be helpful in tuning a PID controlled real-time model 
of beam and ball system, which is an example of double 
integrator system with time delay. 

II. Stability 

Consider the feedback control system shown in Fig. 1. 
The closed-loop transfer function can be written as 
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The characteristic equation of system (2) can be written as 
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By setting <y = cot magnitude and phase of L ( j a ) in 
frequency domain can be written as 
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At gain crossover frequency 

\L ( ja sc) | = 1 
PM =ZL(ja gc ) + 7i 

Substituting (9) into (11) gives 
PM = tan" 1 



By setting K d as a parameter, (20) and (21) can be used to 
choose gain margin (GM) of the system. The stability 
boundary (for gain) can be plotted by setting GM = 1 ■ 
(11) Substituting (3), (4), (5), and (6) into(13), (14), (15), (19), (20), 
and (21) gives 
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Rearranging (12) it can be shown that 

< PM + co er < — 
g 2 

Using (8), (10), and (12) it can be shown that 
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By setting K . as a parameter, (14) and (15) can be used 

to choose phase margin of the system. The stability boundary 
(as far as phase is concerned) can be plotted by setting 

PM = • 

At the phase crossover frequency 

1 
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^L(ja pc ) = n 
Substituting (9) into (17) gives 



(17) 



co pc = tan 



From (18) it can be shown that 
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Using (8), (16), and (18) it can be shown that 
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for phase margin plot and 
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for gain margin plot. 

Combining PM and GM plots for CO that satisfy (22) and 
(25) stability margins for both phase and gain can be observed 

(16) for fixed T , YYl , and K d as on Fig. 2 and 3. Fig. 




(18) 



(19) 
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Fig. 2. Stability margins of the system. Solid line states for gain and 
(20) dashed line for phase. Bold lines are the borders of stability. Below 
both bold lines the system is stable. Each step from higher dashed line 

n/ 

to a lower one increases phase margin by / ■ Each step from 

(21) 

higher solid line to a lower one increases gain margin by 0.2 ■ The 
highest stability margins are obtained in lower left corner. 
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Fig. 3. Stability margins of the system. Solid line states for gain and 
dashed line for phase. Bold lines are the borders of stability. Below 
both bold lines the system is stable. Stability margins are presented as 
on Fig. 2 

2 presents plot for T = 0.35 . K d = 3 and m = 0.7 Fig. 

3 depicts plot for = 0. 1 » K d = 5 and . 

III. Simulation 

Simulation of the control process was made using Matlab 
Simulink software [2-4]. A model was created as on Fig. 1 and 
a simulation was run with 2 sets of parameters. Fig. 4 and 5 
show step and load disturbance responses of the closed- 
loop system for this two sets of parameters. Load distur- 
bance was a 0.2 step. 

1.4 r 




Fig. 4. Step and 
for parameters: 

T = 0.1 • 



20 
time [s] 

load disturbance response of the closed-loop system 

K p = 6.8 , K t = 3.1, K d =5, 



IV. Real-Time Example 



The machine used for real-time experiment was 
Googoltech Ball & Beam digital control system. Transfer 
function of the model of the system can be written as 
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Fig. 5. Step and load disturbance response of the closed-loop system 

K. = 1.3 , K, = 5, 



for 



parameters: 



K 



3.5 



x = 0.1 



T(s) = —re 



(28) 



Control was realized by Matlab Simulink, by an already 
prepared manufacturer model. The user's task was only to 
put PID gain parameters. 
40 r 




time [s] 

Fig. 6. Step response of the real-time system 
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Fig. 7. Step response of the real-time system 



15 



veACEEE 



Poster Paper 



Proc. of Int. Conf. on Advances in Signal Processing and Communication 2013 



Step responses are shown on Fig. 6 and 7, and load 
disturbance response is presented on Fig. 8. 
Experiment was made with parameters: 

K p =6.S (29) 

K i =3.1 (30) 

K d = 5 (31) 
It was observed that results of the experiment differ 
slightly from ones obtained in the simulation. Overshoot in 
the real-time was clearly higher and results were not 
repeatable. The main reason is low quality of the equipment. 
The ball would become stuck on the beam what made the 
model inadequate to reality. Furthermore, the time delay in 
simulation was chosen experimentally and it may be 
inaccurate. 
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Fig. 8. Load disturbance response of the real-time system 
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